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1. Introduction 

I Relativistic effects are important in tfie electrons' dynamics and bound state energies in 

\ lieavy atoms and molecules. When the nucleii involved are heavily charged, the velocities of 

\^ • the electrons of the inner layers are quite large, and so nonrelativistic modelling will lead to 

. important errors. The usual strategies to address this issue are : either use nonrelativistic 

models together with relativistic corrections, or use relativistic models based on the Dirac 

■ operator. 

■ The free Dirac operator is the simplest first order (in time and space) constant coefficient 
, operator which is invariant under the action of the Lorentz group, i.e. its is compatible with 

^ ' Relativity Theory. The free Dirac operator can be written as: 

> : idt+Ha, 
^ ' where 



and the matrices are the Pauli matrices: 



with ai, 


a2 


, Q!3 




' 















1\ /O -A \ 



Two of the main properties of i/c are that its square is a shifted Schrodinger operator: 

-(?h^l\^^m^c^ , (1) 
and that its spectrum is unbounded below and above: 

o-(ifc) = (-00, -mc^] U [mc^, +CX)) , (2) 
and therefore has a gap, (— mc^, mc^). 
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Looking for bound states for an electron in an external potential V of the form 

il>{t,x) = e-'^^ip{x) 

is equivalent to looking for eigenvalues A of the operator + V associated with the 

cigcnfunctions tp. In the case of N electrons, by analogy with the nonrelativistic case, the 

Hamiltonian should be written as : -ffc,xi H 1-, Hc^xn + ^ ■ But as Brown and Ravenhall 

remarked in [5], for potentials V not too singular, the spectrum of this Hamiltonian is 
the whole real line. So, no isolated eigenvalues can correspond to bound state energies. 
But by analogy to the nonrelativistic case, one could still try to find eigenfunctions of 
Hc,xi + • • • +,Hc^xn + V 8iS Slater determinants, that is, in the form; 

V' := ■^^det(V'i(a;j)), i,j = 1,...,N , (3) 

the functions ipj being mutually orthogonal and of norm 1 in the space L2(IR^(C) (the 
constant ^j^,-^ appearing in the above expression just for the sake of normalization). The 

reason to consider Slater determinants is the Pauli principle which states that one cannot 
find two electrons in the same state. Hence, the eigenfunctions should be antisymmetric 
and the simplest antisymmetric function is the determinant. 

The Dirac-Fock equations can be found in the above manner. In that sense, they are the 
relativistic counterpart of the nonrelativistic Hartree-Fock equations which are derived in 
the same way from the Hamiltonian —Ax^ _ . . . _ Ax^^ + V . Despite their non very phys- 
ical derivation, the Dirac-Fock model is widely used in computational atomic physics and 
quantum chemistry to study atoms and molecules involving heavy nucleii. These equations 
were first introduced by Swirles in [29] . One can find many articles in the Physics litterature 
about the Dirac-Fock equations (or more complicated models, like the multi-configuration 
Dirac-Fock model in the case of atoms and molecules with open shells) : see for instance 
[6, 7, 9, 10, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 27, 28]. 

An interesting question is why these equations have bound state solutions (sec below) and 
give very satisfactory numerical results. Why are these equations so well fitted to describe 
the stable electronic configuration in atoms and molecules if the model is not physically 
"well posed" ? In Section 3 we will discuss a possible link between the Dirac-Fock equations 
and models of Quantum field theory which could maybe explain these phenomena. 

A function \1/ e E := = (i/^i, . . . , V'Jv) > (V'^)V'/c)^2 = ^kA is said to satisfy the 
Dirac-Fock equations if for k = 1, . . . N, 

^c^'^k = \ V'fe , Afe e (-1710^,1710^), (4) 

where the mean- field operator H^^ is defined by 

H^^^, = {H, + V)iPk + ip*^)A- [ ^^^^y^^dxdy, (5) 

fI i]R3 \x-y\ 

p being the (scalar) electronic density and the 4x4 complex matrix R, the exchange matrix 
which comes from the antisymmetry of the Slater determinant : 

N N 

Pi^) = ^ {Aix),->Pe{x)) , R{x,y) = ^A{x)<^-tp*{y)- 
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Formally, the solutions of DF equations are the stationary points of the "energy func- 
tional" 

fe=l ^^-^ ^ -^VaxiRS \X-y\ 

in the set 

Remark. The infimum of £c over the set 5* is equal to — oo So, the solutions of DF will 
have to be some kind of saddle points for £c ■ 

2. The Dirac-Fock equations and their nonrelativistic limit 

In an atomic or molecular model, the external potential can be aproximately described by 

V = -Za * , 

\x\ 

where the measure Zjj, represents the nuclear charge distribution, and in the particular case 
of a molecule with K point nucleii with charge Zi and situated at the; points x, .i = 1, . . .K , 
the measure Z^ takes the form Z^i = ZiSx^- So, Z is the total nuclear charge. 

i 

As already said in the Introduction, finding solutions of (4) is then reduced to finding 
critical points of the functional £c on the set S. Once again, the unboundcdness (from 
above and below) of the spectrum of the free Dirac operator makes the functional £c 
totally indefinite. This together with the a priori lack of compactness of the problem ( IR^ 
is unbounded) and the fact that we have to work on a manifold E and not in the whole 
functional space, makes the variational problem difficult. A minimization procedure is once 
again out of the question and another method has to be found. In [11] the authors defined 
a penalized variational problem which could be solved by first maximizing on some part of 
the spinor functions ipi and then defining a more standard min-max argument for a reduced 
functional. 

The theorem proved in [11] states the following : 

Theorem 2.1. ([llj) With the above notations, assume that N and Z are two posi- 
tive integers satisfying max{Z,3N — 1) < ^/2+2/-k '^''^'^ ^ ^ Z + 1. Then, there exists 
an infinite sequence ("^'^'^)^^-^ of critical points of the DF functional £c on E. The func- 
tion vectors (V'l'"', • • ■''P'n) belong to E and they are strong solutions, in H^/'^(JR^ fl 
ni<g<3/2 ^^''^(-''^'^' ^ of the Dirac-Fock equations, that is, for all j > 1, 

H . tpy = el'^ ipy , l<k<N, (6) 



Moreover, 



< e^"'' < ... < e^-'' < mc^ . 

< £c{^'''^) < Nmc^ , 
lim £c{^''^^) = Nmc^ . 



(7) 

(8) 
(9) 
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Remark 2.1. Since is arbitrary, the assumptions of the above theorem contain the case 

of point-like nuclei as well as more realistic nuclear potentials of the form — \^ Pi{x) * i—r , 

F 

where pi G n > , V / Pi = Z . 

Remeirk 2.2. In our units, the above conditions become 

Z < 124 , TV < 41 , N <Z . 

The proof of the above theorem is done by defining a sequence of variational problems 
corresponding to increasing 'topological complexity' or Morse index. By the application of 
our 'first variational argument' wc find solutions 'i''^'^ , that will play an important role (see 
below), since for c large they will be actually 'electronic ground states'. 

Soon after this result was published, E. Paturel used an alternative method to prove the 
same result but without the unnatural condition 

max(Z,3jV-l)< , (10) 

(see [27]). Note that the assumption < Z + 1 is already present in all the existing results 
for the nonrelativistic Hartree-Fock equations. Also, the conditions Z and N less than 
2c^/(27r + 7r/2) are not that unnatural, since already in the linear case such a condition is 
necessary to use Hardy-like inequalities ensuring the existence of a gap of the spectrum of 
Hc + V around (see for instance [30, 31] and [8]). But condition (10) was clearly related 
to the particular method of proof in [11]. Paturel's theorem gets rid of (10) and proves the 
same as Theorem 2.1, but under the sole assumptions 

2c2 2c2 
Z<— ^, N<— -, N<Z + 1. 

7r/2 + 2/7r' 7r/2 + 2/7r' 

Natural questions which arise from the above results are the following: 

- Is there a notion of ground-state in this model? 

- Is there a link between these solutions and the solutions of the usual Hartree-Fock 
solutions? 

- Can we find the solutions in a simple way? By a "simple" variational argument? 

To answer the above questions, it is useful to investigate the nonrelativistic limit of 
the Dirac-Fock equations. Of course, their formal limite arc the well-known Hartree-Fock 
equations, which are identical to the Dirac-Fock ones, but replacing the Dirac operator He 
with the Schrodinger operator — In [12], this was proved rigorously, and as wc sec below, 
this result has been of importance to better understand the variational structure of the 
Dirac-Fock problem and in particular to obtain a good definition of an electronic ground- 
state energy, which is a priori not clear because of the unboundedness of the Dirac-Fock 
energy. 

The following is proved in [12]. 

Theorem 2.2. Let N < Z + 1. Consider a sequence of numbers Cn +oo and a se- 
quence {*"}„ of solutions of (4), i.e. *" = {ip1, • • • , each ip]^ being in H^/^(M!^, (E'^), 
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Assume also that the multipliers £? 



with / V^rV^r^fx = Ski and = e^V^ 

(k = 1, . . . ,N ) satisfy 

< — fij < < ••• < e^v < ~ 'Tiji with fij > nrij > independent of n . 

Then for n large enough, is in H^{W{? and there exists a solution of the Hartree- 
Fock equations, $ = (t^i,--- ,(pn), with negative multipliers, Xi,...,Xn, such that, after 
extraction of a subsequence, 



\l := 4-(c„)^ 



Afc , k = l,...,N , 



(11) 



Xk 







in ^^^(IR^(D^) X ^^^(IR^(D^), 



(12) 



and 



= 0{l/{c^f), 



J:,2(1r3,<c2) 



(13) 



(14) 



Actually the above theorem can be made more precise : the sequence of the solutions 
\]/'=>i given by Theorem 2.1 are shown to converge towards ground state solutions of the 
Hartree-Fock problem when c goes to +oo. This gives a "particular" status to that solution 
of the Dirac-Fock equations. It does not minimize the Dirac-Fock energy among all functions 
in the set S, since inf^ = — oo, but when taking the nonrelativistic limit, they approach 
those solutions of the Hartree-Fock equations which minimize the corresponding Hartree- 
Fock energy. Actually, there is more to it. In [12] the following theorem was proved : 

Theorem 2.3. ([12]) Fix N, Z with N < Z + 1 and take c sufficiently large. Then is 
a solution of the following minimization problem: 



inf {£,(«') ; Gram = , A^ = } 



(15) 



X(-oo,o)(^c,*) is the negative spectral projector of the operator i?c,*> o,nd 



where A,, 



A^f := (A^iPi,--- ,A^xPn) 



Actually, in [12] we call electronic configurations those functions ^ in T, which satisfy 
A^ * = (note that all solutions of the Dirac-Fock equations with positive eigenvalues Sk 
are electronic configurations by definitio). Thus, the above result shows that the Dirac- 
Fock energy is bounded from below in the set of all electronic configurations, and ^"'^ is a 

minimizer for it in that set. 

The above theorems answer some of the questions made above. Another question that we 
asked was that of finding a non very complicated variational argument to find the solutions 
exhibited in Theorem 2.1, since the proofs in that theorem are very involved technically and 
very difficult to implement in actual computations. In order to do this, in [12] we introduced 
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iJc + liJcl I is the positive 



the notion of projector "e-close to A+", where A+ = 
free-energy spectral projector. 

Definition 2.1. : Let be an orthogonal projector in L^(]R^,(D^), whose restriction to 
Hi (m?, (D^) is a bounded operator on {JR^, (C*) . 

Given £ > 0, P+ is said to be e-close to A+ if and only if, for all ip G H's{JR^, <S'^), 



L2(]R.\(i;4) 



< e 



An obvious example of projector e-close to A+ is A+ itself. Other interesting examples 
are the mean-field operators Hc,^ for c large enough. Let us now give a min-max principle 
associated to P+ : 

Theorem 2.4. Fix N,Z with N < Z + 1. Take c > large enough, and P+ a projector 



e-close to A+, for s > small enough. Let P 
E{P+) := inf 



11^2 — P+, and define 



sup 

'I'6(P~//ieSpan(<I>+))" 
Gram^2 



Then E{P+) does not depend on P+ and Eci'^"''^) = E{P+). 

This result shows that when we want to find the electronic ground-state of the Dirac-Fock 

problem by a "simple min-max", we have a large choice of projectors in the limit c — » -|-oo. 
Other characterizations of the ground-state energy for c large can be found in [13]. 



3. The Dirac-Fock equations and a Hartree-Fock model in Fock's space 

Atomic and molecular models based on Q.E.D. seem to be better motivated than the Dirac- 
Fock equations. In a Q.E.D. formulation one should be able to write an energy functional 
which is bounded below in a well chosen subset of the Fock space. In such a model one 
cannot fix the number of particles (electrons) but only the charge, since one has to admit 
the possible presence of electron-positron pairs, which do not change the charge, but change 
the number of particles. In all models based on Q.E.D tehe are some choices to be made, 
for instance, a notion of electron, that is, a positive energy space (and its corresponding 
projector). Also, the kind of configurations that are allowed in the model has to be made 
precise, in order to see if multiconfiguration is allowed or not. 

In [1, 2] such a model has been proposed and the case of no electrons completely treated 

(the empty set). In [4] the same kind of model has been used to treat the case of positive 
charges (A^ > 0). More concretely, in these models one chooses P^ G V, "a positive 
projector" (choice of electronic space), and 

5^.^+ = {7 e 5i(i2) ,7 = 7*, -(1 - P+) < 7 < P+l{l- P+) =0, tr7 < AT} , 

where Si{L'^) is the trace class operators Hilbert space. 

Let us now define 

Kil) = tr {{H. + V - mc^)7) + ^ / ^-^"^^-[^^"J^'''^'' ^xdy , 
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I{c,N,P+)= min jr,(7). (16) 

In a recent result, Barbaroux, Farkas, Helffer and Siedentop have proved that when one 
chooses the class of projectors 

P := {P+, P+ = X[o,oo)(^c,*) , * e E}, 

the minimization problem (16) has a solution. Moreover, imder some (reasonable) condi- 
tions, the minimizer of JF^, in S^'^ is no-pair (purely electronic), that is, there exists 

N 

functions V^i, . . . , V^jv such that : 7^.^ = ^^(V'fe , ■)4'k and 

fc=i 

P+(5,,5)P+^fc = Afc^fe, Afce(0,mc2), (17) 

that is, the 'tp^s are solutions to the projected Dirac-Fock equations (projected in the image 
of the chosen projector P+, of course). 

Now, in [26] , Mittleman proved (somehow formally, but the main ingredients of a proof 

are present in the paper) that any stationary solution (stationary in P+ and 7) of a related 
QED functional would be a solution of the self-consistent Dirac-Fock equations. Hence, 
when we take the energy JF^, to depend both on P"*" and 7 , any stationary point of this 
new functional, (^, P+) , should be a solution of (17) with P+ = x^o +00) i^c,^)- 

A natural idea would be to find stationary points in P+ and 7 by considering the 
max-min problem: 

sup inf jr,(7). (18) 

P+GV 7e5"'-f+ 

This problem was shown to have a unique solution in the case N = in [1, 2]. In the 
case N > , Barbaroux, Farkas, Helffer and Siedentop prove that indeed Mittleman's result 
holds in their context. It is then natural to ask whether the max-min problem defined above 
has a solution or not, and in the affirmative case, whether the solution of that variational 
problem is a stationary point of the energy functional or not. 

In a recent joint work with Barbaroux ([3]), we prove that the answer to the above 
question is yes in some particular cases (for instance when the "neighboring" linear case 
corresponds to a closed shell atom), but it is no in the general case; indeed, we exhibit a case 
where even if the max-min defined above were attained by some (^', P"^) , (^', P^) would 
not be a solution of the self-consistent Dirac-Fock equations, that is P"*" ^ X(o,-i-oo)(-ffc,#) • 

More precisely, it is well known that the Dirac operator is invariant under the action of 
the Lorentz group. This means that for every spatial rotation R, there is a unique matrix 
U{R) e SU2 such that the Dirac Hamiltonian is invariant under the transformation 

and viceversa. 

Then, we prove in [3] the following 
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Proposition 3.1. Assume that fjbc,k (resp. Nc^k) denotes the k-th eigenvalue of the operator 

I 

He + V , not counted with multiplicity (resp. its multiplicity). If N = ^iVc,fe for some 

fe=i 

positive integer I, then for c large enough, the max-min defined in (18) is achieved by a 
solution of the (self-consistent) Dirac-Fock equations, that is, if we denote by P+) a 
solution of (18), then P+ = X(o,+cx3)(-ffc,*) /^'^ ^ large enough. 

Remark. The main reason why the above proposition holds true is that under the as- 
sumptions of the proposition, when c is large enough, for every spatial rotation R, for every 
4* G S, \1/ solution of the Dirac-Fock equations, -ffc,* = -^c, 

I 

In fact we conjecture that when there is no positive integer I such that A'' = Nc^h 

k=l 

(open shell case), the above proposition is not true. Actually, in a special case where V is 

I 

a "small" Coulomb potential, c is large enough and A'' = 1 + Nc,k , I being any positive 

/c=l 

integer, we have proved that Proposition 3.1 does not hold. 

The above results show that there is a clear link between the Dirac-Fock equations and 
some variational problem issued from Q.E.D. since in some cases both methods have the 
same solutions. But this is not the case for every electronic number. 
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